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Abstract: In this paper, we consider a massless field, with spin j, in interaction with 
a Schwarzschild black hole in four dimensions, focusing mainly our study on the s-wave 
scattering. First, using a Fourier analysis, we show that one can have a simple and natural 
description of the Physics near the event horizon without using any conformal field ap- 
proaches. Then, within the same "scattering picture", we derive analytically the imaginary 
part of the highly damped quasinormal complex frequencies and, as a natural consequence 
of our analysis, we show that thermal effects and in particular Hawking radiation, can 
be understood through the scattering of an ingoing s-wave by the non null barrier of the 
Regge- Wheeler potential associated with the Schwarzschild black hole. Finally, with the 
help of the well-known expression of the highly damped quasinormal complex frequencies, 
we propose a heuristic extension of the "tripled Pauli statistics" suggested by Motl, some 
years ago. 
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1 Introduction 

Scattering by black holes (BH) have been studied intensively since the pioneer work of 
Matzner, in 1968 [1]. There have been a lot of works using different approaches. One of 
them uses the complex angular momentum (CAM) theory, the S'-matrix formalism and 
the associated Regge pole techniques. In this framework, the authors have shown [2, 3] 
that the weakly damped quasinormal mode (QNM) complex frequencies of a wide class of 
static, asymptotically flat and spherically symmetric BH of arbitrary dimension, can be 
understood as Breit-Wigner type resonances generated by the interferences and damping 
of a family of "surface waves" lying close to their photon sphere, whose existence is of fun- 
damental importance. It should be noted that, in [4], they have shown that this "surface 
waves" interpretation can also be applied to the BTZ BH, a non-asymptotically flat space- 
time, i.e. in a framework where the notion of an S'-matrix does not exist, by extending a 
powerful formalism introduced several decades ago by Sommerfeld [5]. Finally, by noting 
that each "surface wave" is associated with a Regge pole of the corresponding S'-matrix, 
this approach permits to construct analytically the spectrum of the weakly damped QNM 
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complex frequencies of the corresponding BH, beyond the leading order term. Physically, 
it gives a powerful and quite elegant way to obtain a semiclassical interpretation of BH 
resonance phenomena. However, even if they have computed the absorption cross section 
(beyond the eikonal order) and the greybody factor in such scheme, the WKB approxi- 
mation used to obtain analytically the expression of the Regge poles restricts the study 
to "high frequency" regime, i.e. for frequencies of order w ^ (2Af)~^, where M is the 
mass of the Schwarzschild BH and thus limits the analysis to a scattering description of 
weakly damped QNM, formally associated to high angular momenta. Thus, highly damped 
QNM, Hawking radiation, and more generally, the physics near the BH horizon seems to 
be out of such framework. Nevertheless, in this paper, we claim that one could have a 
physical description of the near horizon physics without using any conformal field theories 
but by focusing, precisely, on the scattering of an ingoing s-wave by the non null barrier 
of the Regge- Wheeler potential of the Schwarzschild BH. In section 2, we introduce some 
notations, assumptions and well-known results linked to a scattering analysis and the S- 
matrix formalism. In section 3, after having discussed the "Rindler approximation", we 
will compute some results which have been assumed or derived by Solodukhin but within 
a conformal field framework [6]. In section 4, we will focus on the "far horizon limit" and 
we will show that one can easily compute the relevant physical quantities needed for a 
scattering analysis. They will allow us to compute the exact expression of the imaginary 
part of highly damped QNM complex frequencies. In section 5, using the results obtained 
in section 4, we will derive explicitly the reflection coefficient related to the scattering of 
the ingoing s-wave by the non null barrier of the Regge- Wheeler potential of the BH which 
is exactly the Boltzmann weight associated with the Hawking radiation, giving, in a purely 
scattering picture, another possible explanation of BH thermal effects. Finally, in section 
6, extending some reflexions due to Motl [7], we will claim by an heuristic approach that 
the real part of these highly damped QNM complex frequencies is intimately linked to some 
"exotic" statistics, going beyond the Bose-Einstein and the Fermi-Dirac ones, which could 
be expected near any BH horizon where a theory of quantum gravity is maybe unavoidable. 



2 Generalities and notations 

We consider first a static spherically symmetric four-dimensional spacetime with metric 

ds^ = -f{r)dt' + ^+r^da^. (2.1) 

Here da^ denotes the line element on the unit sphere S*^. It should be noted that a metric 
such as (2.1) does not describe the most general static spherically symmetric spacetime 
but it will permit us to easily apply the following results to the Schwarzschild case. 
In Eq. (2.1), we shall assume that /(r) is a function of the usual radial Schwarzschild 
coordinate r satisfying the properties: 

• (i) There exists an interval / =]rh, +oo[c R with r/^ > such as /(r) > for r & I. 
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(ii) r/i is a simple root of /(r), i.e. 

f{rh) = and /'(r/,) / 



and 



(2.2) 



These assumptions indicate that the spacetime considered has an single event horizon at 
rh, its exterior corresponding to r € /. If we want to work with an asymptotically flat 
spacetime, in order to define an 5-matrix, we need to impose the condition 



lim /(r) = 1. 



(2.3) 



The Klein-Gordon wave equation for a massless scalar field propagating on a general grav- 
itational background is given by 



1 



-A, 



-g 



0. 



(2.4) 



If the spacetime metric is defined by (2.1), after separation of variables and introduction 
of the radial partial wave functions ^i{r) with £ = 0, 1, 2, ... , this wave equation reduces 
to the Regge- Wheeler equation 



(2.5) 



Here we have assumed a harmonic time dependence exp(— itji) for the massless scalar field. 
The variable r^, = r*(r) is the well-known tortoise coordinate defined, for r € /, by the 
relation dr^/dr = l//(r). Moreover, it is worth noting that the previous definition of r* 
provides a bijection r* = r*(r) from / to ] — oo, -|-oo[. 

In Eq. (2.5), Vi{r) is the Regge- Wheeler potential associated to the massless scalar field. 



V,{r) = f{r) 



(2.6) 



In four dimensions, we recall the well-known cases of non zero spin j fields propagating 
on the BH background described by (2.1). Indeed, one can be interested, for example, in 
the propagation of an electromagnetic test-field (j = 1) or a linearized perturbation of the 
metric {j = 2). Thus, for more general situations in four dimensions, the Regge- Wheeler 
potential can be written as 



Viir) = f{r) 



(2.7) 



where J = 1 — j^, with j the spin of the massless field under consideration. We recall that 
this result is valid when the BH does not carry any electric or magnetic charge. However, 
we would like to refer the reader to [3] for a d-dimensional generalization of the scattering of 
a massless scalar field by a static and spherically symmetric BH through the semi-classical 
CAM techniques. 

According to (2.7), it should be noted that lim.r^rh^ei''') = fiiHr^+oo V£(?") = 
and therefore the solutions of the radial equation (2.5) have a exp(ibia;r*) behavior at the 
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horizon and at infinity. In other words, for a given angular momentum index i, a general 
solution of the Regge- Wheeler equation (2.5), satisfies the following asymptotic behaviors: 



~ (2.8) 

which is a purely ingoing wave at the event horizon and which has the following general 
expression at spatial infinity r* — )• +00 

^u^Ar*) ~ Ane-'^'' + Aoute+'^'' ■ (2.9) 

r»— >+oo 

Moreover, by considering the wronskian of two linearly independent solution of Eq. (2.5) 
at = ±00, we obtain 

l + \Aoutf = \Air,f . (2.10) 

As independent solutions we have chosen Eq. (2.9) and its complex conjugate. Furthermore, 
introducing the transmission and reflection amplitudes and Ri defined by 

T, = ^ (2.11a) 
Ri = ^ (2.11b) 

we can rewrite Eq. (2.10) in the more familiar form 

\Tef + \Re\^ = 1. (2.12) 

Thus, from Eq. (2.9), a QNM, which is defined as a purely ingoing wave at the event 
horizon and a purely outgoing wave at infinity, corresponds to Ain = 0. Therefore, the 
part of the ingoing wave that will not be absorbed by the BH will be reflected back to 
spatial infinity. It should be noted that these quantities can be deduced from the S-matrix 
elements. Here, our problem is spherically symmetric so, the S'-matrix is diagonal. We 
recall that the 5-matrix elements, noted Si{uj), are related to the reflection amplitude by 

Seiu) = i-lY+'Reiu). (2.13) 

It has been shown in [2-4] (and references therein) that the 5-matrix permits to analyze 
the resonant aspects of the considered BH as well as to construct the form factor describing 
the scattering of a monochromatic scalar wave. In this paper, we are mainly interested in 
the scattering of an ingoing s-wave, i.e. i = 0, for which one has, on one hand 

Soiu:) = -i?o(^), (2.14) 

and on the other hand. 

Voir) = ^/'(r)/(r). (2.15) 

For a Schwarzschild BH of mass M, i.e. /(r) = 1 — 2M/r, and a massless scalar field 
(j = 0)) the Regge- Wheeler has a maximum Vo,max for r = (4/3) r/^ [1]. As illustrated 
in Fig. 1, we will show in section 4 that the non null potential barrier, for which we will 
associate an energy Wmin = y^Vo^max) could be at the origin of the Hawking radiation and 
more generally of thermal effects. 
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Figure 1. = 0" Rcgge- Wheeler potential; For the Schwarzschild BH {2M = 1), the maximum of 
the Regge- Wheeler potential is evaluated at Vb,max = (27/256) {2M)-'^ w 0.11 {2M)-'^ and located 
near the event horizon, i.e. at r = (4/3) (2Af ) 



3 The near horizon hmit 



3.1 The Rindler approximation 



In order to deal with the spacetime structure in the very vicinity of the BH horizon, i.e. the 
"near horizon limit", we first naturally expand the function /(r) around the value r = rh- 



f{r)^ f'{rh){r-rh)+ O {r - Vhf . 



Then inserting (3.1) into (2.1), we obtain 

ds^ = -f'f,{r - rh)df + - 
We introduce the new variable p defined by 

dp 

such as 



fh^r - Vh) 



dr 



P 



ir -rh 



or 



r -rh 



Thus, the metric reads 



ds^ 



K^p^dt^ + dp'' + rf^dO. 



2 jr)2 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



which is regular at p = and where k = (l/2)/j^ is the well-known surface gravity of the 
Killing horizon at r = r/j. This "near horizon limit" form of the metric is of course a 
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"Rindler approximation" with a corresponding constant acceleration k. From Eq. (3.1), 
we can obtain the behavior of the tortoise coordinate near the event horizon 

/dr 1 f dr 

which gives 

where the constant of integration is chosen such as r = 2r^ imphes = 0. From (3.7), we 
write the usual radial coordinate r as a function of 

r-rh = rh exp (/^r*) (3.8) 

which allows us to write Eq. (3.1) as 

fin) ~ r,, f'^ exp (/^,) + ^rl /f exp {2f',r,) + ^O^ (exp {3f'^n)) , (3.9) 

where /^^^ = {dP f /drP)\r,^. Then one can express the Regge- Wheeler equation in the "near 
horizon limit" as it has been done in [6], but, as we will see later, it is of no importance 
in this paper. In terms of the tortoise coordinate, the metric (2.1) for the Schwarzschild 
geometry, i.e. /(r) = 1 — 2M/r, in the near horizon limit is given by 

ds"^ = e^^^^* {-dt^ + drl) + r^da^ (3.10) 

where k = 1/4M and r/j = 2M. Of course, if one chooses a hyperplane such as da = 0, 

then in the (t,r^,) coordinates, the metric is conformally flat. 

For the following, we introduce the "null tortoise coordinates" defined by 

u = t + r^, (3.11a) 

v = t-r^ (3.11b) 

and the associated Kruskal coordinates 

kU = e™ (3.12a) 

kV = -e-^". (3.12b) 

3.2 The time dependent "Doppler-gravitational" shift effect 

Now, let us consider, near the Schwarzschild event horizon, a "local inertial frame" with 
coordinates {T,R), i.e. a timelike coordinate T and a spacelike coordinate R, such as the 
Kruskal-Szekeres coordinates are defined as "null coordinates": 

U = T + R (3.13) 
V = T-R. (3.14) 

We recall that the (C/, V) Kruskal coordinate system is analogous to the Rindler coordi- 
nates while the (T, R) coordinate system is analogous to the Minkowski one. Moreover, 
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they cover the entire spacetime manifold of the maximally extended Schwarzschild solution, 
being well-behaved everywhere outside the physical singularity at r = 0. 
Furthermore, since we are working with an s-wave which is spherically symmetric, one could 
always approximate it locally as a "plane wave" propagating in the radial direction. So, in 
the previously defined "local inertial frame" {T,R), we can consider a locally (monochro- 
matic) "plane wave" propagating towards the BH. Up to an amplitude coefficient, it could 
be written as 

^mcrtiai(T, R) a exp [-iLo{T + R)] . (3.15) 

Thus, according to eqs. (3.11) and (3.12), the locally (monochromatic) ingoing plane wave 
^inertial (r,i?) becomes 

$(t, n) = exp \-i (-) e«(*+^')l (3.16) 



in the (t, r*) coordinates where it is obviously not a monochromatic plane wave. It is worth 
noting that keeping the variable u = t + r^, the field reads 



[u ) = exp 



(3.17) 



which is analogous to the "horizon state" (f)H introduced by Solodukhin [6] in the framework 
of a conformal field theory, if one identifies uj/k with the parameter /x/f. Physically, the 
ingoing wave $(n) is not monochromatic simply because of the time dependent Doppler 
shift effect due to the gravitational field of the Schwarzschild BH. 

According to the previous relations, one could conclude that an observer looking at the 
ingoing wave actually sees a superposition of "plane waves" (cf. [8]). Indeed, if we write 
the inverse Fourier transform, we have 

/ + 00 
dOl>o(0)e-*^". (3.18) 
-oo 

The expression of ^{u) can be understood as resulting from a superposition of ingoing 
and outgoing "plane waves" near the Schwarzschild event horizon where <l*o(^^) is the 
amplitude of the "plane wave" of frequency il. In other words, we focused first on a single 
plane wave, i.e. a single frequency. We found that in the freely falling frame it becomes 
a superposition of "plane waves" . Now, if we want to quantize the considered spin j field, 
we have to identify each plane wave with single particle states. In [6], Solodukhin talks 
about a transition between the "horizon state" and the outgoing propagating wave. Here, 
the interpretation from the scattering point of view is analogous and far more natural. 
There is a transition from a single "plane wave" , i.e. one particle state, to a superposition 
of "plane waves", i.e. superposition of single particle states, during the free-fall into the 
Schwarzschild BH, and vice versa by time symmetry. Furthermore, it should be noted that 
close to the horizon, according to eq (2.8), Eq. (3.18) is similar to Eq. (19) in [6] 

/ + 00 
dn$o{n)^n,o{u) (3.19) 
-oo 

but, we can go one step further in this "freely falling wave" description. Indeed, in order 
to know the frequency-distribution of the locally accelerated "wave" ^{u), we can Fourier 
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transform it 



du $(n') e 



/\ JUu' 



(3.20) 



We would like to note that what follows is not really "rigorous" in a quantum field theory 
sense, as it has been noticed in [8], but it permits us to have a physical intuition of what 
happens near the event horizon. According to eqs. (3.17) and (3.20) we can write explicitly 



+ 00 



du' e^^"' exp 



-- 

Thus, following [8], we change variables to y = exp(Kn') and write 

1 r+oo 

^o{n) = - dyy*^/'^-^ exp 
Jo 

The integration reads 



(3.21) 



M^) = -r(—) ( -) '^'^ exp (-7rfi/2K) . 



(3.22) 



(3.23) 



, is naturally 



The modulus squared of the amplitude of each "plane waves" , i.e. 
interpreted as a probability (density) of measuring the associated frequency 0. Since we 
know from the properties of F functions that 



i^l k) sinh (t:VlI K) 



(3.24) 



then the probability (density) of measuring the frequency is given by a Bose-Einstein-like 
distribution 

2 /27r' 



$o(J^) 



1 



(3.25) 



) exp (21:^1 k) — 1 

It is worth noting that this time dependent "Doppler-gravitational" shift effect can be 
adapted to fermions with the formal prescription zO/k— )-zO/k + 1/2. We refer the reader 
to [8] for more details. Of course, the previous quantum point of view with superposition of 
particle states remains valid if one associates the frequency ^ with the energy of a quantum 
particle and interprets the probability of measuring $7 as the probability of finding particles 
of energy 0, as we will show it in section 5. 



4 The far horizon hmit and the QNM complex frequencies 

4.1 The far horizon Hmit 

We start by using the transformation 



1/2 



(4.1) 
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Then, after a partial fraction decomposition, the Regge- Wheeler equation becomes 
1 d\ 



4M2 (ir2 



- + 



■ 2 J -3/4 J + £(£+!)- 1/2 
00^ H H 



r 

8M^u)^ - J - + + 1/2 ^ 4M2a;2 + l/4 

(r - 2M)2 



2Mr 



2M(r - 2M) 



+ 



0. 



(4.2) 



We substitute the variable r by x = {r — 2M)/2M such that x € [0, +oo[ when r G 
[2M, +oo[. Thus, we consider large x region, i.e. r 2M, for which we can expand 
Eq. (4.2) in powers of x~^ 



dx"^ 



+ 



X 



2 ^(^+1) 



1+. 



+ l)-J 



AM'^x^ 



+ 



0. 



(4.3) 



It should be noted that in Eq. (4.3) the spin, J = 1 — j^, only appears in higher orders 
of the expansion, starting from the third one. One way to solve such an equation is to 
consider the first two terms of the asymptotic expansion (4.3) but the price to pay is the 
lost of information concerning the spin 



dx"^ 



+ 



uj^ + + 

X 



0. 



(4.4) 



Thus, the equation becomes analogous to a Coloumb differential equation which reads for 
an s-wave, i.e. £ = 0, 

Vo = 0. (4.5) 



dx"^ 



+ 



a;2 + 



+ 



4.2 The s-vi^ave scattering and the QNM complex frequencies 

In order to find the QNM complex frequencies, we will make some changes of variables. 
We introduce b and z such as 



-4MV = 6(6- 1) 
—AiMojx = z. 

Finally, we define a new "s-wave function" 

il)i_^{z,h) = x'' exp{2iMu;x)(j)uj,oiz). 
Therefore, the equation (4.5) reads 

^dS^ + (25 _ ^)^ _ (5 _ 2iMuj)i;^ 



(4.6a) 
(4.6b) 



(4.7) 



dz^ 



dz 







(4. 



which is one of the well-known confluent hyper geometric equations [9] (for the £ ^ 
case, we refer the reader to [10]). The solution of Eq. (4.8) is a combination of confluent 
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hypergeometric functions of the first kind 



Then the solution for the "s-field" 6^, n reads 



+ (-4iMwx)^-2^ ^b^2iMux - 1) ^(^^i-2b^^ (1 - 6 - 2iML^; 2 - 26; z(x)) 

r(o — 2iMuj) 

(4.10) 

where is defined by Eq. (4.6b). 

For large frequencies, i.e. for a; > Wmin where Wmin has been defined in section 2, we can 
use the following approximation 

b± « l/2±2iMu} 

and, from the large x limit, we use the asymptotic expansions of the previous confluent 
hypergeometric functions in terms of T functions. Moreover, in order to obtain the ampli- 
tudes related to the ingoing wave, i.e. Ain, and to the outgoing wave, i.e. Aout, one has to 
note that x and are related by the identity 

Thus, up to numerical constants of order unity, the amplitude of the asymptotic behavior 
of the s-wave, for "high frequencies" , reads 

-b-2a/^r(l -4fAfw)r(4iMa;) / „i+2b,i-26 



AM . (-4^M.)— ^^^^^^^^^ ( x-n-^^ + 1 ) (4.11a) 



AU-) - (4zMu;)-^^^-^^-^(^-^^"';^^('^^") + l) . (4.11b) 

such as the "s-fleld" has an asymptotical behavior similar to Eq. (2.9). It should be noted 
that, in these expressions, after some simplifications, the value of 6+ = 1/2 + 2iMoj is the 
same for Ain and for Aout- Then, as seen earlier, the QNM complex frequencies which are 
defined by Ain = 0, i.e. by the poles of the function r(l/2 — AiMuj), can be written as 

Vn G N, ^ - 4iMu}n = -n (4.12) 

or, in other words 

Vn G N, STrMuJn = 27ri ( n + - ) . (4.13) 



2^ 

Here, the highly damped QNM complex frequencies are obtained considering a "far region" 
limit, r > 2M, and a "high frequency" regime, i.e. uj > Wmin- It should be noted that, in 
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such calculation, cu^ has no real part, which is known to be non null and depending both on 
the spin j of field and on the characteristics of considered BH. More particularly, one could 
think that not considering the third order term in the expansion in x ^ 1 in Eq. (4.3), is 
equivalent to consider the particular case of spin j satifying 

1- J= l-(l-i') = 0. (4.14) 

In this case, Eq. (4.3) would be equivalent to Eq. (4.5) if j = 0. But we know that this 
"j = result" differs from the results discussed by Motl in [7] obtained with the help of 
the powerful monodromy techniques. In other terms, the statement "Eq. (4.3) would be 
equivalent to Eq. (4.5) if j = 0" seems to be wrong. Then, as we already noticed earlier, 
the exact expression of the "spin j dependent" -QNM frequency spectrum can't obviously 
be obtained by our second order approximation, i.e. Eq. (4.5) for which we only get the 
exact expression of the imaginary part. 

5 From scattering to statistics 

5.1 Quantum field theory: a very brief survey 

In quantum field theory [11], in order to compute the spectrum of outgoing particles from 
a BH, one usually considers the coefficients of the Bogolubov transformation relating the 
well-known creation (resp. annihilation) operator bfi (resp. 6|^) to a^^ and ajj such as 

bn = j du (^a^^a^ - f^Q^al^ (5.1) 
with the commutation relations 

Vio-L,ril2 =5(wi-a;2) (5.2a) 

iVcoi , r/^J = = r]l^, r]l^ (5.2b) 

where r] = a (resp. b) and {001,002) = (to, to') (resp. (0,1^')). The annihilation operators 
6^2 and ajo define respectively the Boulware vacuum state (in the {u, v) coordinate system) 
and the Kruskal vacuum state (in the {U, V) coordinate system) by 

|0_b) = Boulware vacuum (5.3a) 

'^i \^k) = Kruskal vacuum. (5.3b) 

It should be noted that using eq (5.2), the normalization condition for the Bogolubov 
coefficients reads 

I dLo {an.ah,^ - = S {n - n') . (5.4) 

Then, from the standard mode expansions for the considered field operator both in the 
coordinates {u,v) and {U,V), and after some calculations, it follows that the coefficients 
of the Bogolubov transformation obey the well-known relation 

lan^P = exp {SnMn) |/3n^|2 . (5.5) 

Thus, one can easily deduce the expectation value of the "5-particle" number operator, i.e. 
Nq = f'Jj, fen in the Kruskal vacuum state and, more generally, the physics of the Hawking 
effect. 
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5.2 The Hawking radiation: a scattering effect 

Let us focus, now, on the scattering/reflection coefficient of an ingoing s- waves, i.e. Sq{uj) 
—Ro{uj). We recall the main result (4.11) obtained in section 4, i.e. 



AM.) . (4.A,.)-^+^-M.^(l-4^Ma;)^(4^M.) / _,,2.,r-2. + , 

TT V 

Using Stirling's approximation, i.e. considering "large frequencies", i.e. u > comm, one has 
r(l/2 - AiMio) ^ e^iMuiAiMujy^'^^^e-^^^"^ (5.7) 
then, we can easily deduce that the "£ = 0" reflection coefficient can be written as 

i?o(^) = ^ e^iMu. X (5.8) 

From Eq. (5.8), the reflection probability off the = 0" Regge- Wheeler potential barrier 
(cf. Fig. 1) reads 

\Rq{lo)\^ = e-^^^'^^ . (5.9) 

The equation (5.9) could be interpreted as a Boltzmann factor characterized by a temper- 
ature T = I/SttM. With the help of Eq. (5.5), we can now stress the link between the 
quantum field theory approach and the scattering approach, i.e. 

= \Ro{^)\~^ . (5.10) 

5.3 Physics and thermodynamical aspects of the Schwarzschild BH 

In the following, we will associate every "local plane wave" to a "single particle state". 
According to the above results, the infalling motion of a quantum particle is analogous to 
a time dependent boost along the radial direction and consequently the proper distance of 
the infalling quantum particle exponentially decreases with time. Moreover, if the infalling 
particles lose enough energy, they can be considered, for an external observer, as "trapped" 
between the event horizon and the location of Vo^rnax^ which then will play the role of a 
"thermal atmosphere" of the BH. Finally, as a consequence of the time dependent boost, 
this thermal atmosphere will become thinner as the particles eternally fall towards the 
horizon. We claim that this analysis could allow to clarify the link between conformal field 
approaches and BH scattering and could be at the origin of BH thermal effects. Indeed, 
an s-wave, with energy at least equals to i^min, can escape the BH without tunneling. 
Reciprocally, such s-wave will be able to penetrate the barrier from the outside and fall to 
the horizon (cf. Fig. 1). It is worth noting that, in this case, the amplitude of the reflected 
wave off the "^ = 0" Regge- Wheeler potential barrier is exponentially small, but not null. 
Moreover, the mean energy of massless particles in thermal equilibrium at temperature 
T = l/SvrM is roughly of order of T which is bigger than Wmin- Therefore, some of the 
s-waves, i.e. "s-states" particles, will easily escape to infinity. Unless the BH is kept in 
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equilibrium by incoming radiation, it will lose energy to its surroundings. In other words, 
the BH evaporates. This is one possible explanation of the Hawking radiation in terms 
of s-waves scattering, without tunneling. It should be noted that it is not the case for 
fields (or particles) with £ ^ 1 because the Regge- Wheeler potential barrier becomes large 
enough to reduce significantly this process. Moreover, in this case, the location of Vo,max 
moves away from the horizon to be closer to the location of the photon sphere which will 
play a central role in the analysis of weakly damped QNM [3]. 

Furthermore, up to a normalizing contant, the expression (5.9) tells us that the probability 
of finding one particle reflected by the = 0" Regge- Wheeler potential barrier is the 
same as the probability of finding one particle of energy w in a system, i.e. the thermal 
atmosphere, in thermodynamical equilibrium at temperature T = l/SirM. Thus, the 
probability PRiNi^.) to find particles of energy cok in the thermal atmosphere of the BH, 
in thermodynamical equilibrium at the temperaure T = 1 /8ttM is 



1 



(5.11) 



where the normalizing constant Z is the partition function such as the probabilities sum 
up to one 

1 



The mean number of reflected particles depends of course on the nature of the considered 
particles and is usually given by 



1 or oo 

^ Nk=0 

where F/B stands for "Fermion/Bosons" . More explicitly 



(5.13) 



oo 



1 



1 



From the Fourier analysis of section 3, one then has 



for fermions 



for bosons. 



(5.14a) 
(5.14b) 



iV(^/^)(f])\ oc 



(F/B) 



(5.15) 



Therefore, the mean number of particles of energy is proportional to the probability of 
measuring a frequency in the spectrum of the field $(n). 



6 Statistical "heuristics" 

As we have just seen, in thermodynamical terms, the system in thermal equilibrium at 
temperature T = 1 /SvrM is the thin thermal atmosphere located close to the event horizon. 
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In this sense, it is tempting to think that the Physics which has to be considered to 
describe the spectrum of the Hawking radiation would be a conformal field theory, in 
the near vicinity of the horizon, associated to the field in interaction with the considered 
BH. Therefore, due to the accumulation of all the Fourier components (quantum particles 
states) of the incoming wave spectrum on this thermal layer, it would not be a surprise to 
describe the near horizon physics in terms of some "other/new" degrees of freedom which 
would then be expected to have rather "exotic" statistics. In the following, we will give a 
heuristic point of view concerning these new statistics. 

6.1 A naive approach 

Let us introduce the usual variable /3 = to simplify the notations and let us consider 
= 0" fermions reflected off the Regge- Wheeler potential barrier. The mean number of 
fermions with energy a;„ is given by the well-known Fermi-Dirac statistics 

{NiiuJn)) = = (6.1) 

where gi is the degeneracy of the i*^ state of energy uj^- The QNM complex frequencies 
can be seen as poles in the complex cu-plane of the Fermi-Dirac distribution, for which one 
can write 



g/3a;„/T + 1 ~ ^ /3w„ f« 2m + ^ (6.2) 



which is the result obtain in the previous section. 

Even if the behavior is correct for large n, this approach does not give, obviously, the 
correct answer with the well-known constant real part, i.e. ln(3)/87rM. How can one have 
access to the latter? One way to answer this question is to assume that near the horizon, 
there may be another, "exotic" statistics satisfied by the considered quanta which would 
not behave neither like fermions or bosons. 

A first very naive approach is to consider a non-null chemical potential /.i, or a corresponding 
fugacity z = e^^'^ , for a thermodynamical system with a conserved number of fermions. 
Then the Fermi-Dirac distribution reads 

1 z-^ 

(NAoJn - m)) = —3 = —o T. (6.3) 

One could consider the fugacity z as an universal degeneracy g = z~^ for each state of an 
equivalent thermodynamical problem where the number of particles is now not conserved 
and, by definition, which statistics would be given by 

(a,(o.„)) = i (iV,(^„ - ^)) = (6.4) 

Then, the QNM complex frequencies, seen as the poles of the statistical distribution are 

e'^'^" 5 w ^ « In \g\ + 2ni + ^ ■ (6-5) 

The exact analytical result suggest that for a Schwwarzschild BH, 5 = 3. In other words, 
each energy state could have 3 possible sub-levels. We are conscious that such explanation 
is just a mathematical trick and that the underlying physics is still expected. 



-14- 



6.2 More on statistics: the fractional statistics 

However, one could assume that near the event horizon of the Schwarzschild BH, particles 
behave neither like bosons or fermions but they satisfy a more general statistics which 
could be the Polychronakos' fractional statistics [12], the Haldane's fractional statistics 
[13] or even the infinite statistics [14], introduced some years ago, which all start by a 
generalization of the Pauli exclusion principle (highly suggested by the Motl's "tripled 
Pauli statistics" [7]). In the former case, one could use the simple Eq. (22) of [12] which 
reads 

^^^^^-^^)) = ^P^:^- (6.6) 

Fermions and bosons corresponds respectively to a = 1 and a = — 1. If we introduce the 
variable g such as a = 2g + 1, then 

(^-K-g)) = e/^-^ + l + 2g ('-^^ 

where fermions and bosons corresponds now respectively to g' = and g = —1. We deduce 
naturally 

Pujn = \n\l + 2g\+2m(^+^ . (6.8) 

One can noticed that the choice g = 1 (q = 3) allows us to recover the intriguing "tripled 
Pauli statistics" . Let us note that the different expressions of highly damped QNM com- 
plex frequencies can be seen as statistics obtained for different kind of BH which could 
be embedded in the "richness" of these statistics [12-14]. In particular, we think that 
the fractional statistics could be compatible with the "general" expression defining the 
quasinormal frequencies given in [15]. 



7 Conclusion 

In our scattering picture, the Hawking radiation and thermal effects for a Schwarzschild BH 
are associated with the scattering of the s-waves by the "^ = 0" Regge- Wheeler potential 
barrier. A very interesting result is that the scattering/reflection amplitude is exponentially 
small but not null and is exactly the Boltzmann factor related to Hawking effect. We have 
seen, from an external observer perspective, that incoming waves (quantum particles) can 
be "trapped" between the location of Vo,max and the horizon. Moreover, this "thermal 
atmosphere" becomes thinner as a consequence of the time dependent boost, i.e. the free- 
fall of the waves (quantum particles) into the BH. Then one could ask: is there a thickness 
limit for the "thermal atmosphere"? According to the well-known thermodynamical results, 
the answer would probably be linked to the existence of a cutoff, maybe the Planck length, 
but this is way too far from the scope of our analysis. Furthermore, we claim that this 
thin thermal atmosphere could be at the origin of new exotic statistics, due to the change 
of degrees of freedom during the transition between the incoming field in spacetime and 
the accumulation of all its Fourier components (quantum particles states) on the "thermal 
layer". It is worth noting that this thermodynamics, and consequently all the associated 
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quantities, is usually associated to the BH, seen as one thermodynamical system. We have 
thus two different origins for the highly damped QNM complex frequencies. The first is 
the damping of the QNM due to the scattering/reflection off the "£ = 0" Regge- Wheeler 
potential barrier in the "far region limit", directly described by the imaginary part. The 
second would be the thermodynamics of the "thermal layer" associated with the real part. 
Our scattering analysis seems to be in agreement with the Maggiore's interpretation [16]. 
It should be noted that the Hawking effect and the highly damped QNM are not linked to a 
very low frequency limit [1] but simply to the interaction of s-fields with the considered BH, 
whatever their energy is. Moreover, in addition to the existence of a thermal atmosphere, we 
have also recovered some results obtained in [6] but through conformal field frameworks, 
suggesting naturally the use of such approaches to describe the physics near the event 
horizon. Furthermore, one has to bear in mind that the statistics of particles should be 
thought as a property of spacetime and particularly near the horizon of a Schwarzschild 
BH. Therefore, if spacetime is quantum by nature, the commutation or anticommutation 
relations for creation and annihilation operators, for the considered fields, may change 
[17]. Thus, exotic statistics and thereby a modification of the spin statistics theorem could 
be deeply related to the very quantum nature of spacetime and could emerge from an 
underlying theory of quantum gravity in the near vicinity of a BH horizon. Finally, it 
is worth noting that a generalization to d-dimensional Schwarzschild- Tangherlini BH in 
interaction with a scalar (j = 0) field is trivial, at least for the imaginary part of the QNM 
complex frequencies. Indeed, in this case, the "dimension parameter" in the Regge- Wheeler 
potential is multiplied by the angular momentum £ (for more details, see [3]). Thus, for the 
scattering of an s-wave, the results will be identical. Concerning the real part, there is no 
"simple" argument except that the "£ = 0" Regge- Wheeler potential is not modified and 
thus one could perform the same calculations and recover the same scattering/reflection 
amplitude as in section 4. Actually, the real part seems to not be affected either by the 
dimension of spacetime as it has been noticed by Motl in [7] (and references therein). 
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